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Abstract 

Efforts to give an improved mathematical meaning to Feynman's 
path integral formulation of quantum mechanics started soon after 
its introduction and continue to this day. In the present paper, one 
common thread of development is followed over many years, with con- 
tributions made by various authors. The present version of this line 
of development involves a continuous-time regularization for a gen- 
eral phase space path integral and provides, in the author's opinion 
at least, perhaps the optimal formulation of the path integral. 

The Feynman Path Integral, 1948 

Much has already been written about Feynman path integrals, and, no doubt, 
much more will be written in the future. A comprehensive survey after 
more than fifty years since their introduction would be a major undertaking, 
and this paper is not such a survey. Rather, it is an attempt to follow 
one relatively narrow development regarding a special form of regularization 
used in the definition of path integrals. Since we deal with several different 
approaches, this paper does not go too deeply into any one of them; it is 
intended more as a conceptual overview rather than a detailed exposition. 
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In order to set the stage, let us start our discussion with a review of the 
traditional approach to path integral construction. 
We begin with the Schrodinger equation 

..dib(x,t) H 2 d 2 ib(x,t) < \ i / \ /-x 
at 2m ox 2 

appropriate to a particle of mass m moving in a potential V(x), x e M.. A 
solution to this equation can be written as an integral, 

i>(x",t") = J K(x",t";x',t')4>(x',t')dx' , (2) 

which represents the wave function ip(x", t") at time t" as a linear superposi- 
tion over the wave function ip(x', t') at the initial time t', t' < t" . The integral 
kernel K(x", t"; x', t') is known as the propagator, and according to Feynman 
[1] it may be given by 

K(x", t"; x', t')=N J e «/R)/[(W2)* a (*)-v(x(t))] * Vx ? (3) 

which is a formal expression symbolizing an integral over a suitable set of 
paths. This integral is supposed to run over all continuous paths x(t), t' < 
t < t", where x(t") = x" and x(t') = x' are fixed end points for all paths. 
Note that the integrand involves the classical Lagrangian for the system. 

Unfortunately, although highly suggestive, the preceding expression for 
the path integral is undefined as it stands: For example, the normalization 
constant Af diverges, and the putative translation invariant measure Vx does 
not exist. To overcome these basic problems, Feynman adopted a lattice 
regularization as a procedure to yield well-defined integrals which was then 
followed by a limit as the lattice spacing goes to zero called the continuum 
limit. With e > denoting the lattice spacing, the details regarding the 
lattice regularization procedure are given by 

K(x", t"; x', t') = \im{m/2Tiihtf N+1 ^ 2 J - J 

x exp{(i//i)l£ [(m/2e)(x l+1 - Xl f - eV(x t )]} IT^ dx, , (4) 

where x N+1 = x", x = x', and e = (*" - t')/(N + 1), N G {1,2,3, . . .}. In 
this version, at least, we have an expression that has a reasonable chance of 
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being well denned, provided, of course, that one interprets the conditionally 
convergent integrals involved in an appropriate manner. One common and 
fully acceptable interpretation adds a convergence factor to the exponent of 
the preceding integral in the form 



which is a term that formally makes no contribution to the final result in 
the continuum limit save for ensuring that the integrals involved are now 
rendered absolutely convergent. 

Accepting the fact that the integrals all converge ensures that a meaning- 
ful function of e emerges, but, by itself, that fact does not ensure convergence 
as e — > and, even if convergence holds, it furthermore does not guarantee 
that the result is correct! To ensure convergence to the correct result requires 
that some technical condition(s) must be imposed on the potential V(x). In 
this regard, we only observe that a correct result emerges whenever the po- 
tential has a lower bound, i.e., whenever V(x) > c, for some c, — oo < c < oo, 
for all x. 

We recall that for the free particle of mass m the potential V(x) = for 
all x. In that case, Eq. (4) reads 



which is the form of the quantum mechanical propagator for the free particle. 
Comments 

The procedure sketched above — whereby the action functional expressed 
as an integral over a continuous-time parameter is replaced by a natural 
Riemann sum approximation, which eventually is followed by a continuum 
limit (e — > 0) as the final step — provides a satisfactory procedure for a 
suitable and large class of potentials to define the propagator K(x", t"; x', t'). 
However, it is important to stress that a lattice regularization followed by a 
continuum limit is by no means the only way to give a satisfactory definition 




(5) 




x exp{(i/h)^ =0 [(m/2e)(x l+1 - Xl ) 2 } U^ =1 dx l 
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of a path integral, nor, in the author's opinion does it even represent the 
most satisfactory definition, although admittedly this latter issue is in part 
subjective. What follows in this paper is a narrow review — an historical 
slice — of the development of various efforts to find a suitable continuous- 
time regularization procedure along with a subsequent limit to remove that 
regularization that ultimately should yield the correct propagator. While 
some of the work to be described did not follow directly from work that 
preceded it, all of this work does, with appropriate hindsight, seem to have 
a fairly natural progression that makes an interesting story. 



Feynman-Kac Formula, 1951 

Through his own research, Mark Kac was fully aware of Wiener's theory of 
Brownian motion and the associated diffusion equation that describes the 
corresponding distribution function. Therefore, it is not surprising that he 
was well prepared to give a path integral expression in the sense of Feynman 
for an equation similar to the time-dependent Schrodinger equation save for 
a rotation of the time variable by —tt/2 in the complex plane, namely, by the 
change t — > —it. In particular, Kac [2] considered the equation 

^ = ^-V(. W ,,«). (7) 

This equation is analogous to Schrodinger's equation but of course differs 
from it in certain details. Besides certain constants which are different, and 
the change t — > —it, the nature of the dependent variable function p(x, t) is 
quite different from the normal quantum mechanical wave function. For one 
thing, if the function p is initially real it will remain real as time proceeds. 
Less obvious is the fact that if p(x, t) > for all x at some time t, then 
the function will continue to be nonnegative for all time t. Thus we can 
interpret p(x, t) more like a probability density; in fact in the special case 
that V(x) = 0, then p(x,t) is the probability density for a Brownian particle 
which underlies the Wiener measure. In this regard, v is called the diffusion 
constant. 

The fundamental solution of Eq. (7) with V(x) = is readily given as 



4 



which describes the solution to the diffusion equation subject to the initial 
condition 

\im + W(x,T;y,0) = 5(x-y) . (9) 

Moreover, it follows that the solution of the diffusion equation for a general 
initial condition is given by 

*VO = /wV,t>(*',0^. do) 

Iteration of this equation N times, with e = (t" — t')/(N + 1), leads to the 
equation 

p(x", t") = N' j ■■■ J e -W^(Lo(xi + i-*i) 2 Uf =1 dxi p(d, t') dx' , (11) 

where and xq = x' . This equation features the imaginary time 

propagator for a free particle of unit mass as given by 

W(x", t"; x', t') = N' J ■■■ J e -(V2^)s^ =0 (, i+1 -, i ) 2 u^dxi . (12) 

Since this equation holds for all iV > 0, we may assume that it also holds in 
the limit e — > 0, i.e., N — > oo, which we can write formally as 

W(x", t"; x', t')=N J e -^I ±2 dt Vx , (13) 

where M denotes a formal normalization factor. (Symbols such as M may 
stand for different factors in different expressions.) 

The similarity of this expression with the Feynman path integral [for 
V(x) = 0] is clear, but there is a profound difference between these equa- 
tions. In the former (Feynman) case the underlying measure is only finitely 
additive, while in the latter (Wiener) case the continuum limit actually de- 
fines a genuine measure, i.e., a countably additive measure on paths, which is 
a version of the justly famous Wiener measure. In particular, 

W(x",t";x',t') = J d^ w (x), (14) 



where ^ denotes a measure on continuous paths x(t), t' <t < t", for which 
x(t") = x" and x(t') = x' . Such a measure is said to be a pinned Wiener 
measure, since it specifies its path values at two time points, i.e., at t — t' 
and at t — t" > t'. (The traditional Wiener measure, which we shall not 
really deal with in this paper, specifies its values only at the initial time, and 
it corresponds to one additional integration over the final value x" at the 
final time t" .) 

We note without proof that Brownian motion paths have the property 
that with probability one they are concentrated on continuous paths. How- 
ever, it is also true that the time derivative of a Brownian path is almost 
nowhere defined, which means that, with probability one, x(t) = ±00 for all 
t. 

When the potential V(x) ^ the propagator associated with (7) is for- 
mally given by 




an expression which is well defined if V(x) > c, — 00 < c < 00. A mathemat- 
ically improved expression for (15) makes use of the Wiener measure and is 
given by 

W(x", t"; x', t') = J e -fvW))* d^ w {x) . (16) 

This is an elegant relation in that it represents a solution to the differential 
equation (7) in the form of an integral over Brownian motion paths suitably 
weighted by the potential V. Incidentally, since the propagator is evidently 
a strictly positive function, it follows that the solution of the differential 
equation (7) is nonnegative for all time t provided it is nonnegative for any 
particular time value. 

Gel'fand and Yaglom, 1956 

In an important review article [3], these two authors introduced the concept 
of a continuous-time regularization of Feynman path integrals. Although 
their procedure was later shown to be incorrect, their work may be said to 
have initiated an interesting line of development. 
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The idea of Gel'fand and Yaglom was relatively straightforward. Formally 
stated, their proposal was to define the propagator as 



lim M v [ exp{(i/h) f[(m/2)x 2 - V(x)]dt} exp{-(l/2i/) fx 2 alt} Vx . (17) 

i/-»oo J 

Formally, therefore, their proposal consisted of introducing an auxiliary fac- 
tor into the integrand that is identical to the factor which leads to Wiener 
measure. The purpose of such a factor is to introduce a regularization into 
the original formal expression. Finally, the limit v — > oo is taken as a last 
step which amounts to formally removing the auxiliary factor and leaving 
behind the original integrand of the Feynman path integral. 

There is one glaring difficulty with this proposal. The auxiliary factor 
which leads to Brownian motion paths with a diffusion constant v < oo 
has the property of ensuring that the paths are continuous, but they are 
nowhere different iable. For such paths, the contribution of the potential is 
well defined, but the contribution of the kinetic energy is divergent for all 
paths. Gel'fand and Yaglom hoped to get around this difficulty by combining 
the kinetic energy with the corresponding factor in the auxiliary term giving 
rise to a Wiener measure with a complex diffusion coefficient a, where 

- = --^. (18) 
a v n 

Thus the strategy was to take the limit of a sequence of presumably well 
defined integrals over Wiener measure with a complex diffusion constant a 
as the real part of a vanished. Notice that this kind of regularization does 
not involve the introduction of a temporal lattice followed by a continuum 
limit but rather maintains a continuous time parameter throughout. It is 
for this reason that we refer to this kind of procedure as a continuous-time 
regularization scheme. 



Cameron, 1960 

Shortly after the appearance of the Gel'fand and Yaglom paper, Cameron [4] 
showed that the scheme was fatally flawed. In particular, Cameron showed 
that a Wiener measure defined with a complex diffusion constant a leads to 
a countably additive measure only when a = Re<r > 0, hence Imcr = 0. It 
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is instructive to sketch the argument that leads to this conclusion. Consider 
the iV-fold integral 

/ A \ (w+i)/2 r r 

{£-) J- J exp[-(A/2e) (x l+1 - xtf] dx t 

= (^y-) exp[-(A/2iV e ) (x N+1 - x ) 2 } , (19) 

which represents one of the primary formulas involved in discussing Wiener 
measures. This formula holds for any complex A so long as Re A > 0. If 
such a formula exists in the limit that e — > and iV — > 00 such that iVe 
remains constant, then the desired path measure will also exist. Although 
that condition would appear to be evident, it is not true for general A. 

Integrals exist when they converge absolutely. Absolute convergence cer- 
tainly holds for (19) for all N < 00, but what about when N — > 00? To test 
this issue, we consider the absolute integral associated with (19) given by 

/ IAI \ ( N + 1 )/ 2 r r 

(ht) J "J ex PH ReA / 2e > s £o(^+i - *tf] n,1, dx, 

= (j^H^f exp[-(ReA/2iV £ )(, N+1 - . (20) 

Evidently, the final result exists as iV — > oo only provided |A| = Re A > 0, 
i.e., provided ImA = 0, as claimed. Consequently, the proposal of Gel'fand 
and Yaglom fails on the ground that whenever the diffusion constant for the 
Wiener measure is truly complex (Imcr 7^ 0), then the so-defined Wiener 
measure is only finitely additive and not countably additive. In this case, 
finite additivity means that finite answers for general integrals involving this 
measure arise by delicate cancellation of both positive and negative divergent 
contributions. A regularization scheme that results in a finitely additive mea- 
sure on paths is not qualitatively different from the original formal Feynman 
path integral, and therefore it does not constitute an acceptable continuous 
time regularization. 

ltd, 1962 

Soon after Cameron pointed out the lack of a countably additive measure 
for the Gel'fand- Yaglom procedure, Ito [5] proposed another version of a 
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continuous-time regularization that resolved some of the troublesome issues. 
In essence, the proposal of Ito takes the form given by 

lim M v J ex^{{i/K)j[\mx 2 - V{x)\dt) exp{-(l/2u) f[x 2 + x 2 ]dt} Vx .(21) 

Note well the alternative form of the auxiliary factor introduced as a regu- 
lator. The additional term x 2 , the square of the second derivative of x, acts 
to smooth out the paths sufficiently well so that in the case of (21) both 
x(t) and x(t) are continuous functions, leaving x(t) as the term which does 
not exist. However, since only x and x appear in the rest of the integrand, 
the indicated path integral can be well defined; this is already a positive 
contribution all by itself! 

To proceed further with (21) we need to decide on what values to fix at 
the initial and final times. As was the case with the procedure proposed 
by Gel'fand and Yaglom, let us first suppose that a composition law of the 
expression in (21) holds before the limit v — > oo is taken. This requires 
that one fix not only x(t") = x" and x{t') = x' but, in addition, that we 
fix x(t") = x" and x(t') = x'. Viewed as a solution of the Schrodinger 
equation, however, we are at a loss to understand the meaning of the values 
of x and x to be held at any time. Thus, while this interpretation may yield 
a limiting function it is hard to see that such a result could be a solution to 
Schrodinger's equation. 

In point of fact, Ito did not choose the previous set of data to be fixed but 
he chose another set. Ito interpreted (21) as an integral in which x(t") = x" 
and x(t') = x' are the only values held fixed. He therefore did not require 
any composition law on the expression (21) at finite v, but only expected a 
composition law to hold for the expression that arises after v — > oo. 

With such an interpretation in mind we proceed to evaluate in a fairly 
heuristic manner the functional integral (21). For simplicity we shall discuss 
only the simple case where V(x) = 0; moreover, we set t' = 0, t" = T, and 
choose x(0) = 0, x(T) = x. Thus, prior to taking the limit v — > oo, let us 
focus on the formal integral 

Af u J exp{( y i/h)fxg( y t)dt}exp{-(l/2u)f[x 2 + a 2 x 2 ]dt}5(x(0))Vx, (22) 

where a 2 = 1 — imu/h, g(t) will be chosen below, and for the foregoing 
expression we assume there is no other pinning; the required second pinning 
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at t" = T will be introduced in a moment. Let us set 

*(t)s (23) 
for all t, so that the former integral can be replaced by 

Kj exp{(i/^)/^(*) dt} exp{-(l/2u)J[i 2 + a 2 ?) dt}V^ . (24) 

In this form, the integral involves an Ornstein-Uhlenbeck measure [6], and 
the answer is readily given by 

exp[-{u/Aah 2 )jg{t)g{u)e- alt - ul dtdu] , (25) 

where we have chosen M v so that the answer is unity if g(t) = for all t. As 
a next step we let g(t) = A for < t < T, and g(t) = otherwise. Thus (25) 
becomes 

exp[-(z/A 2 /4a^ 2 )/o/o e " a|t " M| dtdu \ ■ ( 26 ) 

Since 



J e W»JWt)* e -W»* dX /(2irh) 



it follows that 



5(x(T)-x), (27) 



A^y 5(x(r) -x) e -(wJ\i 2 +- 2 e]dt v ^ 

= J e -W Xx e - v>?F ' iah2 d\/(2nh) 

a ( ax 2 \ 

vFi ^ ( " Tf) ■ < 28) 



where 



2T 2 
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As v becomes large, we may set a 2 ~ —imu/h, in which case a/vF may be 
replaced by —im/2Th. Hence, we are led to the final result 

/ m fimx 2 \ 

VM exp (w)' (30) 

which is recognized as the quantum mechanical propagator K(x, T; 0, 0) for 
a free particle of mass m as given in (6). 

Ito [5] made this story rigorous for constant, linear, and quadratic poten- 
tials, as well as those potentials of the form 

V(x) = Je txs w(s)ds (31) 

provided that f\w(s)\ds < oo. 

In summary, by introducing a regularization involving a higher derivative 
(i; 2 ), Ito was able to soften the paths sufficiently to allow the kinetic energy to 
be well defined. In so doing he was able to give a satisfactory continuous-time 
regularization for (21) for a certain class of potentials V. 



Feynman, 1951 

It is necessary to retrace history at this point to recall the introduction of 
the phase space path integral by Feynman [7]. In Appendix B to this article, 
Feynman introduced a formal expression for the configuration or g-space 
propagator given by 

K(q",t";q',t')=M J exp{(i/ft) f[pq - H(p, q)]dt} VpVq . (32) 

In this equation one is instructed to integrate over all paths q(t), t' < t < 
t", with q(t") = q" and q(t') = q' held fixed, as well as to integrate over 
all paths p(t), t' < t < t", without restriction. As customary, this is a 
formal statement and in practice it needs to be given a precise meaning. 
The lattice prescription in which the continuous time parameter is replaced 
by a finite set of discrete points is the procedure that is typically followed. 
For completeness, we illustrate a common lattice space version of the formal 
phase space path integral expression, as given by 

K(q",t";q',t') = hm j ••• J exp{(i//i)E^ [|p m/2 (^ +1 - q t ) 

-eH( Pl+ y 2 ^(q l+1 + q l ))}}Ul d Pl+1/2 /(27rh)UfL 1 dq l . (33) 
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In this expression, all p and q variables are integrated over except for qN+i = 
q" and go = q', and, just as before, e = (t" — t')/ (N + 1). Since q\ implies a 
sharp q value at time t' + le, we have chosen to name the conjugate variable 
Pi+1/2 to emphasize that a sharp p value must occur at a different time, here 
at £' + (/ + l/2)e, since it is not possible to have sharp p and q values at 
the same time. Note that there is one more p integration than q integration 
in this formulation. This discrepancy becomes clear when one imposes the 
composition law which requires that 

K(q>", f- q\ = J K{q'\ t">; q", t") K{q\ t"; q\ t') dq" , (34) 

a relation which implies, just on dimensional grounds, that there must be one 
more p integration than q integration in the definition of each K expression. 

Observe that (32) exhibits an apparent covariance under canonical coor- 
dinate transformations, but this is quite illusory. In fact, the lattice form 
(33) is correct only when the canonical variables are described by Cartesian 
coordinates, and this kind of limitation on straightforward lattice-space reg- 
ularizations is quite general. In a later section of this paper we shall return 
to this point and explain why this limitation is necessary. 

It is also instructive to consider an alternative phase space path integral 
for the momentum or p-space propagator formally given by 

K{p\t"-p\t')=M j exp{(i/h)f[-qp-H(p,q)]dt}VpVq, (35) 

which is obtained from (32) by Fourier transformation on both end variables, 
q" and q'. In this case a lattice space definition can be given by 

K(p",t";p',t') =lim J ••• J e W {(i/mlo[-h+iMPi+i-Pi) 

-eH{\{p l+1 +p l ),q l+1/2 )\} nf =1 d Pl Uf =0 dq l+1/2 /(2nh) , (36) 

and we see in this expression that there is one more q integration than p 
integration. Again this makes sense when one considers the composition law 

K(p>", t'"; p', = J K(p>", r ; p", t") K(p", t"; p', f ) dp" . (37) 
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On the other hand, and following similar reasoning, one would be hard 
pressed to give a satisfactory lattice space formulation of the putative formal 
phase space path integral given by 

M J exp{(i/h)f[±(pq-qp)-H(p,q)]dt}VpVq. (38) 

(In the light of the previous discussion, even the physical meaning of such an 
expression seems uncertain, at least to the present author. On the other hand, 
such an expression is easily understood in the coherent state representations 
that follow.) 

It is widely appreciated that the phase space path integral is more gen- 
erally applicable than the original, Lagrangian, version of the path integral. 
For instance, the original configuration space path integral is satisfactory for 
Lagrangians of the general form 

L(x) = \mx 2 + A{x)x -V(x) , (39) 

but it is unsuitable, for example, for the case of a relativistic particle with 
the Lagrangian 

L( x ) = -mVl - x 2 (40) 

expressed in units where the speed of light is unity. For such a system — as 
well as many more general expressions — the phase space form of the path 
integral is to be preferred. In particular, for the relativistic free particle, the 
phase space path integral 

interpreted in the sense of (33), is readily evaluated and yields the correct 
propagator. 

Issues of proper coordinate choice also arise for the configuration space 
path integral as well, and expressions such as (3) and (4) implicitly refer to 
Cartesian coordinates. 

Coherent State Representations, 1960 

As a prelude to the following section it is pedagogically useful to recall some 
basic properties of coherent states and the Hilbert space representations they 
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generate [8]. We focus on coherent states generated by Heisenberg canonical 
operators P and Q which obey the basic commutation relation [Q, P] = ihl, 
where I denotes the unit operator. More specifically, we shall assume the 
Weyl form of the commutation relation holds for self-adjoint momentum and 
position operators in the form 

e -iqP/h e ipQ/h _ e -ipq/2h e i(pQ-qP)/h = jj^^ ^ 

where p and q are arbitrary real variables. In an abstract Hilbert space, let 
us introduce a normalized fiducial vector \rj), which is otherwise arbitrary, 
and define vectors of the form 

\p,q) = U[p,q]\r ] ) (43) 

for all (p,q) G R 2 . These states are the canonical coherent states, and they 
have a number of interesting properties. Since U[p, q] is a unitary operator, 
it follows that each vector is normalized, i.e., 

||b,?>ll = <P,?b,9> 1/2 = l- (44) 

Next, it should be noted that the coherent states are continuously parame- 
terized; specifically, if (p, q) — > (p 1 , q') in the sense that \p— p'\ 2 + \q — q'\ 2 — > 0, 
then it follows for any choice of \rj), that || \p, q) — \p', q') || — > 0, or as one says, 
the vectors are strongly continuous in the parameters p and q. Finally, we 
want to stress the important property that these vectors admit a rather con- 
ventional looking resolution of unity as a superposition over one-dimensional 
projection operators onto the coherent states, and specifically that 

J \p,q)(p,q\dpdq/(2nh) = I (45) 

holds, where, as before, / is the unit operator, and this relation holds for any 
choice of \rj). It is easy to show that the resolution of unity holds weakly in 
the sense that, for arbitrary vectors \ip) and \<f>), 

J {<j>\p,q){p,q\^)dpdq/(2nh) = . (46) 

(It is less easy to show, but nevertheless true, that the resolution of unity 
also holds as a strong operator identity.) 
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(42) 



The resolution of unity formula permits us to introduce novel representa- 
tions of Hilbert space rather different from those customarily used, particu- 
larly in quantum mechanics. Let us introduce the functional representatives 

ijj(p, q) = (p, q\ip) , (47) 

which are defined for all in the Hilbert space. As we have already done for 
the coherent states themselves, we suppress any dependence of the functions 
ip(p,q) on the state Indeed, for any choice of each such function is 
a bounded, continuous function of the variables p and q. For each pair of 
such functions, such as ip(p,q) = (p, q\i[>) and (fi(p,q) = (p,q\(j>), we assign 
the inner product given by 

(0, V) = J <f>(p, g)>(p, q) dpdq/(27rh) = {m . (48) 

The last part of this relation shows that the functional inner product equals 
the inner product in the abstract Hilbert space. The so-defined functional 
representation is complete simply because every element of the abstract 
Hilbert space is imaged as a function of the form ip(p,q). It should be no- 
ticed that the space of functions defined in this way is a complete and proper 
subspace of the space of square integrable functions over two variables. The 
space of functions {ip(p,q) = (p, q\ip)} defined for all \ip) and all (p, q) con- 
stitutes the coherent state functional representation of a one-particle Hilbert 
space. In this representation, for example, the basic Heisenberg operators 
are given by 

P ~ iH dq ' Q^q + ^— , (49) 

and thus the coherent state representation of the Schrodinger equation is 
given by 

ihdt/j(p : q, t)/dt — H(—ihd/dq, q + ihd/dp)tp(p, q, t) . (50) 

What makes this a one-particle problem (rather than a special kind of two- 
particle problem) is the proper choice of acceptable initial conditions for this 
equation. In particular, it suffices to take as an initial condition the coherent 
state overlap function (p, q\p', q'), for arbitrary p' and q'. 
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The solution to this form of Schrodinger's equation can be expressed in 
the form 

*P(p",q"X) = J K(p",q",t";p',q',t')ii(p',q',t')dp'dq'/(2nh) , (51) 

where K(p", q", t"; p', q', t') denotes the propagator in the coherent state rep- 
resentation. 

The coherent state representation differs from the usual Schrodinger rep- 
resentation given in (1), and it also has a different physical interpretation as 
well. For one thing, the meaning of ^(a;)) 2 is the probability density to find 
the particle at position x. The meaning of \ip(p, q)\ 2 , however, is rather dif- 
ferent; this quantity is simply the probability that the state \ip) can be found 
in the state \p, q). The physical meaning of the variable x in the Schrodinger 
representation is that of a sharp position. Correspondingly, we may ask what 
is the physical meaning of the parameters p and g? To answer that question 
it is convenient to restrict the choice of the fiducial vector very slightly by 
imposing the conditions that {77 1 1 77) = {r)\Q\rj) = 0. In that case, it follows 
that 

{p,q\P\p,q) =p , {p,q\Q\p,q} = q , (52) 

implying that unlike x in the Schrodinger representation, the variables p and 
q represent mean values in the coherent states. It is for this reason that 
we can specify both values simultaneously at the same time, something that 
could not be done if instead they both had represented sharp eigenvalues. 

The propagator for the coherent state representation of Schrodinger's 
equation can also be given a formal phase space path integral form, namely 

K(p",q",t";p',q',t')=M J exp{(i//i)/[pg - H(p,q)]dt} VpVq , (53) 

which superficially is just the same expression as (32)! What makes these 
expressions different is what values are pinned and more explicitly what are 
the respective lattice space formulations. In the coherent state case, we have 

K(p",q",t";P,q,t') = lim J ••• J exp{(i/^)Ejl [i(p J+1 + Pl )(q l+1 - qi) 

-eH(l(pi +1 + pi) + i\(qi+i - q t ), \{qi+i + qi) - i\{pi+i - Pi))]} 
x exp{-(l/4ft)££ [(p m -p{f + (q l+1 - gi ) 2 ]} Ul 1 dp l dq l /(27fh)(5A) 
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Observe that there are the same number of p and q integrations in this 
expression. Such a conclusion is fully in accord with the combination law as 
expressed in the coherent state representation, namely 

K{pi\q>\t'"-p>,q>,t>) 

= j K(p>", q">, f- p\ q\ t") K{p\ q\ t"; p> , q>, f ) dp"dq"/(2nh) .(55) 

Daubechies and Klauder, 1985 

By any measure, phase space is the natural arena for both classical and 
quantum mechanics. In classical physics, dynamics evolves in phase space 
in a highly symmetric fashion for general Hamiltonians, while in quantum 
physics both P and Q appear in remarkably similar roles within the general 
formalism. Moreover, as stressed previously, the phase space path integral is 
more widely applicable than the original configuration space path integral. 

In 1985, Daubechies and Klauder [9] examined the phase space path in- 
tegral in a new light. In particular, they were led to consider the expression 

Jim^y* exp{if\pq-H(p,q)]dt} exp{-±f[p 2 + q 2 ]dt}VpVq, (56) 

which extends the continuous-time regularization procedure to phase space 
path integrals. Let us initially look at this expression from a general view- 
point. Evidently we can interpret the regularization as two independent 
Brownian motion regularizations in the limit as the diffusion constant di- 
verges, which is rather like the original Gel'fand-Yaglom procedure. Al- 
ternatively, and since from a classical point of view p is often rather like q 
(depending on H, of course), this procedure has some features superficially in 
common with that of Ito. However, there are also important differences with 
those previous procedures as well. The chosen regularization ensures that 
both p and q are continuous functions but that p and q are almost nowhere 
defined. The derivatives enter the integrand only in the form jpqdt. How- 
ever, this term can be interpreted as f pdq, which for Brownian motion paths 
is a standard stochastic integral that has been well studied in both its Ito 
and Stratonovich versions; for reasons to be given below, when it becomes 
important to choose a rule we shall choose the Stratonovich rule. We inter- 
pret the integration in (56) to be pinned for both p and q at both end points: 
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specifically, (pit") , q(t")) = {p",q") and ip(t') , q(t')) = (p',q')- In this case it 
is clear that the result of (56) yields a function of the form 

K(p",q",t";p',q',t'). (57) 

By the time of 1985, representations of propagators for one-particle quan- 
tum systems in the form of coherent state representations, such as those de- 
scribed in an earlier section, were well known to many researchers. If, indeed, 
(56) turned out to yield coherent state representations for solutions to the 
Schrodinger equation — as opposed, say, to either the g-space or p-space rep- 
resentations — there would be no reason for concern. In fact, that is exactly 
what happens! 

In [9] it was rigorously established that 

lim M v [ *tM--ff(P.5)]* e -( 1 /2^j "\p 2 +q 2 ]dt J) p J) q 
v^oo J 

= lim 2nhe uT/2h [ e {i/h) $ [pdq ~ H{p ' q)dt] d^ w (p,q) 

v^oo J 

= (p",q"\e-^"-^\p>,q>) 

= K(p",q",t";p>,q',t'), (58) 

where the second line of (58) is a mathematically rigorous formulation of 
the heuristic and formal first line, and /i^ denotes pinned Wiener measure. 
In fact, the connection indicated here is far reaching in that the choice of 



regularization also determines that 

\p,q) =e- iqP/h e ipQ/h \0) , (59) 

[Q,P]=ihI, (Q + iP)|0) = 0, (60) 

H = J H(p,q)\p,q)(p,q\dpdq/(27ch) . (61) 

A sufficient set of technical assumptions ensuring the validity of this repre- 
sentation is given by 

a) J H(p, qf e ~ a{p2+q2) dpdq < oo, for all a > , (62) 

b) J H(p, qfe-^ p2+q2) dpdq < oo, for some (3 < l/2h , (63) 

c) 7i is essentially self adjoint on the span of finitely many 
number eigenstates . (64) 
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A wide class of Hamiltonians — but by no means all acceptable such 
operators — consists of all semibounded, symmetric (Hermitian) polynomials 
of P and Q. We also note that the connection between the operator H and 
its symbol H(p, q) given in (61) is that of anti-normal ordering. In particular, 
this means that 

e -(r-is)(Q+iP)/2H e (r+is)(Q-iP)/2H = J e i( Sg -rp)/h^ ^ ^ g | d pdq/{2lth) , (65) 

and any other connection follows by suitable linear combinations. As an 
alternative association we also note that 

H(p, q) = e -(«/4)(<M) H w {p, q) , (66) 

where Hyy(p, q) denotes the well-known Weyl symbol of the operator Ti. [In 
most of the author's earlier work, the particular symbol H(p, q) used here 
has been denoted by the lower case symbol h(p, q).] 

Observe that the auxiliary factor in (56) that provides the regularization 
involves a metric, specifically, da 2 = da(p, q) 2 = dp 2 + dq 2 , on classical 
phase space. Such a metric characterizes a flat, two-dimensional phase space 
expressed in Cartesian coordinates. 

It is interesting to consider the transformation of (58) under a canonical 
change of phase space coordinates. To this end we introduce p = p(p, q) and 
q = q(p, q) as two new classical canonical coordinates that are related to the 
original canonical coordinates by the relation 

pdq = pdq + dF(q, q) , (67) 

where F is known as the generator of the transformation. This is a standard 
form for this kind of relation in which q and q are regarded as the independent 
variables; another version that we will find more useful is given by 

pdq + dG(p, q) = pdq . (68) 

Under such a canonical coordinate transformation, the metric da 2 = dp 2 +dq 2 
assumes the form 

da 2 = da(p, q) 2 = A(p, q) dp 2 + 2B(p, q) dpdq + C(p, q) dq 2 , (69) 
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for suitable A, B, and C, but, nevertheless, the metric still characterizes 
the same flat, two-dimensional phase space even though it is now expressed, 
generally speaking, in curvilinear canonical coordinates. 

In the original Cartesian coordinates the stochastic integral fpdq eval- 
uated by the Stratonovich rule is identical to its evaluation using the It 6 
rule. Under coordinate transformations, we adopt the Stratonovich defini- 
tion [mid-point rule, cf.(54)] since for this rule the differential and integral 
laws of ordinary calculus still apply. Consequently, the differential rules (67) 
and (68) which apply for classical functions (say C 2 ), apply equally well 
to Brownian paths (which are only C°). Recognizing that the Hamiltonian 
function transforms as a scalar, we see that 

H(p, q) = H(p(p, q), q(p, q)) = H(p, q) . (70) 

Assembling the separate parts, we learn that under a canonical coordinate 
transformation, (58) becomes 

lim M u /" e (V^/[p^(p,9)-^(p,9)]^ e -(i/2^)/[da(p,g) 2 M 2 ]^ x>pVq 

= lim 2Trhe uT/2h [ e m)J\P<z^m-Hm<*] dji u w (p,q) 

= @\q"\e- i{t "- t,)n/h \p',q') 

= K(p",q",t"-p',q',t>), (71) 

where ~]i v w denotes pinned Wiener measure on a two-dimensional flat phase 
space expressed in curvilinear coordinates. The term 

JdG(p, q) = G(p", q") ~ Gtf, q') , (72) 

and thus its presence amounts to no more than a phase change of the corre- 
sponding coherent states. Specifically, 

\p^q) = e - iG (p^)/ h e - i P(p^)Q/ h e i i(p^) p / h ^ — e~ iG ^'^^ h \p,q} , (73) 
[Q,P]=ihI, (Q + iP)|0) = 0, (74) 

H = I H(p,q) \p,q)(P,Q\ dpdq/(2nh) . (75) 

Under a canonical coordinate transformation, therefore, and apart from a 
possible phase change, observe that the coherent state vectors are unchanged; 
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it is only their labels that have changed. Observe further that the Hamiltonian 
operator H is completely unchanged; only the details of its representation 
are different in the new coordinate system. 

In summary, we see that a continuous-time, Brownian motion regulariza- 
tion of the phase space path integral can be rigorously established. It applies 
to a wide class of Hamiltonians, and the formulation is fully covariant under 
general canonical coordinate transformations. The only additional ingredi- 
ent to attain this covariance is the introduction of a metric, da 2 , on classical 
phase space which is ultimately used to support the Wiener measure regu- 
larization. As a general rule, it is noteworthy that quantization schemes that 
agree with orthodox quantum mechanics — and therefore agree with a vast 
number of experiments — do not exhibit manifest covariance under canoni- 
cal coordinate transformations. Hence, the fact that the present procedure 
exhibits covariance under canonical coordinate transformations is especially 
to be welcomed. 

As a final comment, we remind the reader of the general requirement that 
canonical quantization, i.e., the choice of certain phase space variables to be 
"promoted" to operators, should be carried out in Cartesian coordinates 
[10]. This rule has caused much concern over the years and its elimination 
has been, partially at least, a motivating factor in certain alternative quanti- 
zation schemes. Instead of seeking to eliminate the Cartesian coordinates, we 
fully accept them, and, indeed, we understand the need for them as follows 
[11]: With only a symplectic form on phase space — and in the absence of 
any phase space metric — there is no self-referential mechanism to associate 
the correct physical meaning of any given mathematical expression, say, for 
the Hamiltonian H(p, q). In particular, if H(p,q) = p 2 /2, and in the absence 
of any further information, how is one to determine that this particular ex- 
pression is the Hamiltonian for a free particle, or for an harmonic oscillator, 
or for some particular anharmonic oscillator, etc., any one of which can be 
brought to that mathematical form by a suitable canonical coordinate trans- 
formation. It is with the addition of a flat-space phase-space metric, and with 
due regard to its coordinatized expression, that the physical significance of a 
given mathematical expression can be determined, and with that determina- 
tion, the quantization procedure can lead to the correct energy spectrum for 
the physical system under consideration. The phase space continuous-time 
regularization scheme discussed in this section offers the great advantage that 
the introduction of the metric into the formal phase space path integral si- 
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multaneously gives mathematical and physical meaning to the path integral 
by assuring that the proper interpretation of the Hamiltonian is maintained 
throughout the quantization procedure. In point of fact, one may say that 
the procedure we have illustrated offers a true geometric quantization proce- 
dure, and consequently, the entire procedure can therefore be expressed in a 
coordinate free form. That is indeed the case, and the interested reader is 
referred to [12] for further discussion of this point. 

Final remarks 

A number of further developments and applications regarding continuous- 
time regularizations have appeared in recent years. We do not enter into any 
details of these works here, but rather list just a few additional papers for 
further study: [13, 14, 15]. 

Dedication 

It is a distinct pleasure to dedicate this paper to my long time friend and 
colleague, Hiroshi Ezawa. I have known Hiroshi for approximately forty years 
and I have profited greatly from the many times we have been together and 
during which we have often worked closely on a variety of topics of mutual 
interest. I wish him many more years of good health and productive research. 
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